Magnetic Properties of HTSC with Weak Interlayer Coupling. 
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Magnetic properties of layered high temperature superconductors with a weak interlayer coupling in 
the region of critical fluctuations are investigated in framework of the Ginzburg-Landau approach. 
The sample magnetization is calculated perturbatively to the second order in the interlayer coupling 
constant. The first order correction can be incorporated into the non-interacting expression for the 
magnetization with the corresponding shift of the critical temperature. Only the second order con- 
tribution modifies thermodynamics, leading to violation of the 2D scaling law and to disappearance 
of the magnetization crossing point. Magnetic field suppresses the interlayer interaction, and at 
sufficiently high values of the applied field 2D scaling is recovered. 



I. INTRODUCTION. 

Scaling properties of layered high temperature super- 
conductors (HTSC) around the mean-filed transition line 
H C 2(T) has Japen under intensive experimentaloa and 
theoreticalHily study during the last decade. The lavr 
ered YBa2Cu307_^ compounds exhibit 3D properties!. 
All thermodynamic and transport quantities obey the 
3D scaling law as a function of the scaling variable 
B(T-T c2 (H))/(TH) 2 / 3 , when the applied magnetic field 
is perpendicular to the layers plane. n In the case of 
Bi2Sr2Ca2Cu30io experimental resultsB show 2D scal- 
ing behaviour with respect to the 2D scaling variable 
A(T — T c2 (H)) / (TH) 1 / 2 . The fact that the phenomenon 
is observed in the vicinity of the mean-field upper criti- 
cal field indicates that the problem can be studied in the 
framework of the Ginzburg-Landau (GL) theory with the 
order parameter projected onto the space of the lowest 
Landau levels (LLL). First calculations of ths-,£see en- 
ergy scaling function were done perturbativelyETQ. The 
non-perturbative approach |br pure 2D systems was de- 
veloped by Tesanovic et. al.u. It was shown that around 
the mean-field transition line only the fluctuations of 
the total amplitude of the order parameter play impor- 
tant role in the superconductor thermodynamics. The 
remaining part, fluctuations in the position of the vor- 
tices enters through the dimensionless AbrikosovEH geo- 
metric factor /3a- In the most realistic cases this quan- 
tity weakly depends on the vortex configuration and can 
be taken as a constant from the very beginning. This 
approach was generalized for thadayered 2D and 3D sys- 
tems by Tesanovic and AndreevEI. They showed that the 
system dimensionality depends on the strength of the 
interlayer coupling constant. When coupling is absent, 
the system becomes pure two dimensional. In the limit 
of strong interlayer coupling, when the superconducting 
correlation length £ c along the c-axis direction becomes 
much larger then the effective interlayer separation s the 
3D description becomes appropriate. In the intermediate 
case £ c < s the system becomes quasi two-dimensional 
for which scaling is impossible. 

In this Paper we study quasi two dimensional layered 
superconductor in the region of the critical fluctuations. 



In this region we use GL-LLL description, which allows 
us to build a solvable model. The interlayer coupling is 
assumed to be small and therefore can be taken into ac- 
count perturbatively. Within this approach we calculated 
the sample magnetization to the second order in the in- 
terlayer coupling constant. If this constant is extremely 
small one can consider only the first order correction. Wc 
show that because of the nearest neighboring coupling 
this correction leads to shift of the critical temperature, 
preserving the form of usual 2D scaling function of the 
magnetization. Only the next order correction destroys 
scaling and leads to disappearance of the magnetization 
crossing point. The effective interlayer coupling constant 
turns out to be proportional to the square root of the 
inverse applied field, which results in effective suppres- 
sion of the interlayer interaction in high field region and, 
therefore, leads to recovering of the 2D scaling law. 



II. THE MODEL. 

Consider a layered type II superconductor in the region 
of critical fluctuations around its mean-field transition 
line T c2 (H) (or H c2 {T)). The applied magnetic field H 
is assumed to be normal to the layers plane: H || c. Then 
the superconductor thermodynamics at the temperature 
T can be described by the following partition function: 



Z oc J J P[A] expj- 



•T 7 !*, A] 
knT 



(1) 



where J 7 ^, A] is the GL free energy functional of the lay- 
ered system with a nearest neighboring Josephson cou- 
pling between pancake vortices given by 



J 7 !*, A] = sJ2 [ d 2 r |ao|* 

12 , |,t, ,T, 12 , (H- B„) 2 



2 + fl*n| 4 - 
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(2) 



and fee is the Boltzmann constant. The quantity s is 
an effective interlayer spacing, ao — a(T — T c q) and 
(3 = const are the first and the second GL coefficients 
correspondingly. The third GL coefficient 7 is assumed 
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to be anisotropic, where the quantities j a b and 7 C define 
its value in the layer plane and c direction correspond- 
ingly. In what follows, we refer to the quantity j c as 
an interlayer coupling constant. The quantity 9 n (r) is 
the order parameter of the nth layer and B„(r) || c is the 
magnetic induction induced in the nth layer of the super- 
conductor. Two-dimensional gauge invariant gradient d n 
is defined as: 



-it) 



d 2e 
dr c 



A„(r), 



where A„(r) = V x B„(r). The sample magnetization is 
given by 
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where Nl is the number of layers. In the weak coupling 
regime, it is convenient to introduce renormalized critical 
temperature: 

T c = T c0 - 2^ = T c0 - ST C . (4) 
a 

This modifies the expression (^) for the free energy, which 
now reads as 



lab\d n 9 n \ 2 ~ 7c*£(#n-l + *n+l) 



| 2 + fl*n| 4 - 

(H-B ra ) 2 
8tt 



, (5) 



where a = a(T — T c ). Thus, in order to solve the prob- 
lem one has to calculate extremely complicated integrals 
over the order parameter 'J and the vector potential A 
appearing in the expressions (^) and (^|) for the parti- 
tion function and the magnetization. However, there are 
number of simplifications, which can be done. 

In the limit of large values of the GL parameter k 
(ft ~ 100 for the most HTSC) one-pan neglect fluctua- 
tions of the magnetic induction B„E3. Then we minimize 
the last expression (^) for the GL free energy with respect 
to the vector potential A n as it is done in the case of con- 
ventional superconductor. This leads to the set of decou- 
pled (with respect to the layer index) GL equations for 
the vector potential A„ and the order parameter *& n (r). 
In general, the order parameter can be expanded over 
the electron eigenfunctions of the Landau levels in the 
applied magnetic field H. However, close to the mean- 
field transition line H C 2{T) = —a(T)c/(2hej ac ) one can 
restrict oneself to the zeroth Landau level only. This ap- 
proximation can be used, at least, if H > 1/3H C 2(T), 
until the first Landau-level becomes important. Re- 
cently, it was shownli§E-3 that the lowest Landau level 
(LLL) approximation works good even if H <C H C 2(T) 



for k ^> 1. In the LLL approximation the equations, 
described above, can be solved analyticallytil. After sub- 
stituting these solutions into the expression (0) for the 
free energy, we finally obtain: 



T = sS 



a 1 



H 

Hc2 



9, 



7 C ** (*„_!+ * 



n+l) j 



(6) 



where the bar means averaging over the layer area S. In 
this case only the integrals over the order parameter * 
remain in the expression for the partition function (uh: 



Z (x V[9] exp 



k R T 



where the quantity J-^] is given now by (||). In the 
LLL approximation the sample magnetization (|^) can be 
calculated using the following formula: 
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87raft 2 L ^ 



J V[9] exp - 



(7) 



In order to proceed further, we define the Abrikosov 
geometric factor for each layer separately: (3a (n) = 

W^/(W^)\ and its average f3 A = A^E/^n). 

Following the refsBEl we assume that the quantity /3a (n) 
only slightly depends on the actual vortex configuration 
and therefore we put /3a(1) = /3a (2) = •■• = Pa- Then 

we replace the quantities \9 n \ A by /3a (j^nl 2 ) in the ex- 
pression (^|) for the free energy. With this replacement 
the model becomes exactly solvable: 



Z oc / T>[A] exp < -N v 



x\A r 



4 (' A "' 2 ) : 



- M A* (A„+i + A„_ x )] } , (8) 



where N v — ^/^o is the number of vortices. The stan- 
dard 2D scaling variable x is given by 



x = A 



T-T c2 {H) 



where A = y/s^o/(16TTK 2 p A kB)H' c2 and H' c2 = 
—dH C 2(T)/dT\rr—T c . The dimensionless coupling con- 
stant u is 



u = A- 



5T C 



2VTH 

and the dimensionless order parameter A„ reads as 

2/3a/3 



(9) 



|A n r = A- 
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With these new variables the sample magnetization (]?]) 
can be expressed as 



M 



a AT— i 



A 



(10) 



III. 



CALCULATION OF THE PARTITION 
FUNCTION. 



In order to compute magnetization of the superconduc- 
tor we, first, have to calculate the partition function (||). 
This involves evaluation of the integrals over the order 
parameter A. The main difficulty in such calculation is 
the quartic term appearing in the exponent of the right 
hand side of the formula (||). In order to decouple it, we 
introduce a set of additional integration variables {A„}: 

Z oc J V[A] J V[\] exp [ A « + 

(x + iA„)|A„| 2 - /xA* (A„+i + A„_!)] } , (11) 



where 



r„ 

The contours T n are parallel to the real axis, standing on 
some distance from it, in order to insure convergence of 
the integrals over the order parameter A. The formula 
( |ll| ) is the result of use of the simplified version of the 
Hubbard-Stratonovich transformation, usually applied in 
the field theory. As it was explained above, the order 
parameter A n is the linear combination of the electron 
eigenfunctions of the lowest Landau level: 



A «( r ) ^^C nk L k {v), 



fe=0 



and 



exp 



2/ 2 



where I is the magnetic length corresponding to the 
charge 2e. Then the meaning of the integration over the 
order parameter becomes clear: 



V[A] cx J]/ dC* nk dC nk . 

n,k 



The integrals over these expansion coefficients in (|ll| ) arc 
of the generalized gaussian type and can be evaluated 
analytically. As a result, we obtain: 



Z(x J V[X]exp{-N v C} . 



(12) 



The action C is given by 



C = tr 



A 2 + In ^xl + iX — /ij^j 



where A is diagonal matrix, consisting from the elements 
A„. The matrix 7 is the real symmetric matrix, arising 
as a result of the interlayer coupling: 

The integrals over the expansion coefficients C nk con- 
verge, if along the integration contours r„ the following 
inequality is satisfied: 



Re(/„) > 0, 



(14) 



where /„ are eigenvalues of the complex symmetric ma- 
trix xl + iX — /17. In the thermodynamic limit N v —* 00 
the integrals over A„ in the partition function (|l^) can 
be evaluated within saddle point approximation. It will 
be shown below that due to particular properties of the 
saddle-point manifold the condition (|l4|) can be satisfied 
in the weak-coupling regime in rather large range of the 
scaling variable x. 

The saddle point equation is found from the condition 
SC = and reads as 



tr 



2A + i I xl + iX — jj/) 



SX > = 0. 



(15) 



The position of the saddle point depends now on the tem- 
perature T and the applied field H via two parameters 
x(H, T) and u(H, T). In this case, 2D scaling of the mag- 
netizationrjjlg) becomes impossible, as it was predicted 
in the rem. If the coupling is weak 2\l < \x + iX\ Xthis 
inequality is similar to Tesanovic-Andreev criteriorD for 
quasi 2D systems), the second term in the left hand side 
of the equation ( Jig ) can be expanded in the powers of /i. 
Then the saddle point equation (|l5| ) can be rewritten as 
follows: 



2X r 



X + iX n 
2x + i(Xn-l + A„+i) 



(x + iX n ) 2 (x + iX n -l)(x + iXn+l) 



+ o(fi 4 ) = 0. (16) 



The terms of the order of [i and /i 3 drop out in the right 
hand side of the last equation, since tr(7 2n+1 ) = 0. The 
structure of the saddle point equation is such that, at 
least, up to the second order in /i the saddle point solution 
for A is proportional to the unit matrix, namely 



Ai — A2 



A. 



Then the equation ( |T6| ) can be rewritten in the following 
simple form: 



2A 



2ia 2 
1 

iX (x + iX) 3 



+ oGu 4 ) = 0. 



(17) 



3 



As soon as the saddle point solution is found, the sample 
magnetization ([!(]) can be calculated using the following 
simple relation: 



Aki 



M 



2Ak B 



iX. 



(18) 



From the last equation we conclude that the physically 
meaningful solution for A lies on the imaginary axis. Fur- 
ther, using the Rayleigh-Ritz theoremEj, one can show 
that both matrices 21 — 7 and 21 + 7 are positive defi- 
nite. Then the eigenvalues of the matrix 7 belong to the 
interval [—2,2]. In this case, the conditions (|l4|) for con- 
vergence of the integrals over expansion coefficients C^ n 
in the partition function can be written as: 



2[i < \x + iX\, 



(19) 



which is automatically satisfied in the weak coupling 
regime. 

IV. MAGNETIC PROPERTIES. 

In the previous section we derived the saddle point 
equation (|17]) for the layered superconductor under as- 
sumption that the interlayer coupling constant fi is small. 
The sample magnetization M(H,T) is proportional to 
the saddle point solution and is given by the formula 
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FIG. 1. The quantity Mj\fWf as a function of the 2D 
scaling variable x is plotted for 8T C — IK. 

In the " zeroth" approximation only one of the two sad- 
dle point solutions of the eq. (|l7|) can be reached by 
allowed deformation of the integration contour: 



X (x) = -(x- y/x 2 +2). 

Then the magnetization preserves the 2D scaling law and 
is given by the same formula as in the noninteracting 
caseO: 



M 



— (x— \fi 



2). 



(20) 



Actually, this "zeroth" order result includes the first or- 
der correction in /j, by means of the critical temperature 
shift ([|) . The second order correction to the saddle point 
solution is calculated as a small perturbation: 



X(x,n) = Xo{x) 1 + 4/i 



2 «A (x) 



Vx 1 



Then the sample magnetization is given by 



M = M 1 - 2/x 2 



V ' x 1 



(21) 



(22) 



where Mq(H,T) is the magnetization of the "decoupled" 
sample given by the equation (pc|). Using the last for- 
mula, we plotted the quantity M/ V HT as a function of 
2D scaling variable x (see the figure [l]) for five differ- 
ent values of the applied field H between 10 and 50kOe. 
The interlayer coupling constant is chosen to be small, 
such that 5T C = lK(see the eq. || for definition of ST C ). 
We used T c = 111K, n = 100, H' c2 = iOkOeK" 1 and 
s = 2nm, which aee of the order of the typical exper- 



imental parametersa. 
variable x in the fig 
dition of the theory (|l9) 



The whole range of the scaling 



satisfies the applicability con- 
. As it was expected, the 2D 
scaling is destroyed, due to the dimcnsionless coupling 
constant /i(T, H) appearing in the right hand side of the 
expression (|22j) for the magnetization. 
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FIG. 2. The sample magnetization M as a function of the 
temperature T is plotted for 5T C = IK. 

The plot of the magnetization as a function of tempera- 
ture can be found in the fig. g, in which we used the same 
set of the phenomenological parameters, as for the fig. [|. 
Like in the previous case, the whole temperature range in 
this figure satisfies the applicability condition ([l^) . Form 
this figure follows that even account of the small inter- 
layer coupling destroys the magnetization crossing point, 
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which is the property of the scaling form (see ref.Ea) . 
Actually, the crossing point is noticeably destroyed only 
in the low-field region. The higher the applied field, the 
better the crossing point is pronounced. This is the con- 
sequence of the specific form of the effective coupling con- 
stant /x (see eq. The interaction correction to the 
magnetization in the formula ( p^ ) is proportional to the 
square of this constant and therefore is proportional to 
H^ 1 . Then, at the large values of the applied field in- 
terlayer interaction is effectively suppressed. Indeed, one 
can treat the magnetization data plotted in the figure ^| as 
if they were obtained from the experiment. Then, we try 
to fit them to 2D scaling form variating phenomenological 
parameters T c and H' c2 , as it is usually done in experi- 
ments. The fitting results are given in figure [| The beat 
fit is obtained for T c « 111.5K. Like in the experiments, 
the scaling results are insensitive to the value of the phe- 
nomenological parameter ^oH' c2 in relatively large range 
of its values around ATK . It can be observed from 
the figure || that the scaling is satisfactory good for fields 
larger than 30kOe and is destroyed in the weak-field re- 
gion. It must be stressed that the mentioned above fit is 
made essentially by hand and without any error estima- 
tion. However, it can serve as demonstration of effective 
suppression of interlayer interaction in the high field re- 
gion. 



to reduce the problem to the simpler one, namely to the 
saddle point equation (0). This equation can be solved 
perturbatively. We calculated the magnetization of the 
superconducting sample to the second order in the effec- 
tive coupling constant /i(T, H). If the coupling is suffi- 
ciently weak, one can consider the first order correction 
only. It turns out that this correction does not modify 
the scaling properties of the sample, leading to the trivial 
renormalization of the critical temperature (see eq. (|])). 
Account of the second order correction leads to violation 
of the 2D scaling and destroys the magnetization crossing 
point in the low field region. At sufficiently high values 
of the applied field the interlayer interaction is effectively 
suppressed. This leads to recovering of 2D scaling with 
a well pronounced crossing point. 
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FIG. 3. 2D scaling of the magnetization data, taken from 
the fig. |. 



V. SUMMARY. 

In this paper we considered HTSC assuming weak in- 
terlayer nearest neighboring coupling. In the region of 
the critical fluctuations, close to the mean-field transi- 
tion line H C 2(T) the order parameter can be taken as a 
linear combination of the electron eigenfunctions of the 
lowest Landau level. This approximation together with 
the assumption that the Abrikosov geometric factor only 
weakly depends on the actual vortex configuration allows 
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